Asymmetric electron-phonon interactions in the three-band Peierls-Hubbard model 
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Using the Quantum Monte Carlo (QMC) technique within frozen-phonon, we studied the ef- 
fects of the half-breathing O(7r,0) phonon mode on the ground-state properties of the three-band 
Peierls-Hubbard model. Our simulations are performed for both ionic and covalent electron-phonon 
couplings. The effects of lattice displacements on the ground-state energies and charge fluctuations 
are similar in magnitude for both hole- and electron-doped cases. However, the effects of lattice 
displacements on the magnetic properties are rather different. In the hole-doped case, the nor- 
malized next-nearest-neighbor Cu-Cu spin correlations are dramatically modified by both ionic and 
covalent electron-phonon couplings. On the other hand, in the electron-doped case, much smaller 
effects are observed. The distinct spin-phonon couplings, in conjunction with the spin-bag picture 
of the quasiparticle, could explain a strong mass renormalization effect in the p-type cuprates and 
a weaker effect in the n-type cuprates. 
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Recently, a large number of experimental observa- 
tions in the superconducting cuprates suggest that the 
electron-phonon (el-ph) coupling in these systems is 
strong and -in conjunction with strong electronic corre- 
lations, may play a role for the superconductivity mech- 
anism. In particular, angle resolved photoemission spec- 
troscopy (ARPES) data indicate a mass renormalization 
( "kink" ) in the electronic dispersion of the p-typc high-T c 
cuprates near a characteristic energy scale (ss 70 meV)i£, 
which is possibly caused by coupling of quasiparticle to 
some phonon modes. Moreover, ARPES data 3 also show 
that this mass renormalization is weak in the n-type 
Ndi.85Ceo.i5Cu04, indicating that the el-ph coupling in 
the p-type materials is much stronger than in the n-type 
ones. On the other hand, both inelastic neutron scat- 
terin g 4 i 5 i 6 and inelastic x-ray scattering^ found that the 
anomalous softening of the Cu-0 bond-stretching oxygen 
mode is present in both p- and n-type materials, suggest- 
ing that although there exist some differences, the el-ph 
coupling is strong in both kinds of materials. 



In order to understand this asymmetry of el-ph 
coupling in p- and n-type cuprates, we have studied 
the three-band Peierls-Hubbard model in the physi- 
cally interesting parameter regime by employing the 
constrained-path Monte Carlo (CPMC) methodi2*±i. In 
the CPMC method, the ground-state wave function \^q) 
is projected from a known trial wave function \^>t) and 
then the physical quantities are calculated in this ground 
state. We find that in both hole- and electron-doped 
cases, the effects of lattice displacements on the ground- 
state energies and charge fluctuations are similar. How- 
ever, an essential difference between hole- and electron- 
doped systems is that the response of spin degree of free- 
dom to phonons is strong in the former case, but rather 
weak for the latter. 



In the hole representation the three-band Peierls- 



Hubbard model 8,9 has the Hamiltonian, 
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Here the operator d ia creates a hole at a Cu 3d x 2_ y 2 or- 
bital and p'jp creates a hole in an O 2p x or 2p y orbital. 
Ud and U p are the Coulomb energies at the Cu and O 
sites, respectively. V p d denotes the nearest-neighbor (nn) 
Coulomb repulsion. For the el-ph coupling, we consider 
that the nn Cu-0 hybridization is modified by the O- 
ion displacement Uj as t p J d = <pij(t p d ± ctUj), where the 
+ (— ) applies if the bond shrinks (stretches) with pos- 
itive Uj. The phase factor <^>y takes a minus sign for 
j = i + x/2 and j = i — y/2. The Cu-site energy is as- 
sumed to be modulated by the nn O-ion displacements 
Uj as e d ({Uj}) = €d + PJ2j^ zU j- The other electronic 
matrix elements are the 0-0 hybridization ij£ = ±i pp 
with the minus sign occurring for k = j — x/2 — y/2 and 
k = j +x/2 + y/2, and the O-site energy e p . The charge- 
transfer energy is defined as e = e p — td- For simplicity, 
we introduce dimensionless parameters: u = a\uj\/t p d, 
or u — f3\uj\/t p d- In units of t p d, we use the set of pa- 
rameters: Ud = 6, e = 3, U p = Vpd — 0, and t pp = 0.5. 

All the results reported here were done on a lattice 
with 6x6 unit cells and periodic boundary conditions. In 
the following we will focus on the half-breathing O(7r,0) 
phonon mode and examine the dependence of physical 
quantities on the lattice displacement u in different cases. 
In order to diminish systematic and statistical errors in 
the CPMC simulations, we studied two closed-shell cases: 
the hole doping Sh — for the hole-doped case and 
the electron doping S e = |g for the electron-doped case. 
Notice that we were forced to choose two different dop- 
ings for the hole- and electron-doped cases, due to the 
requirement of a closed-shell condition. However, we 
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FIG. 1: Displacement pattern of the half-breathing O(7r,0) 
mode on the Cu-0 plane. The closed circles label the Cu 
sites, and the open circles the O sites. The arrows denote the 
lattice displacements from the O sites. 



have checked that the results at a smaller electron doping 
S e — on a lattice with 8x8 unit cells remain qualita- 
tively similar to the ones at 5 e = demonstrating that 
our findings on the 6x6 lattice reflect intrinsic properties 
of the three-band Pcicrls-Hubbard model. 



Fig. shows the displacement pattern of the half- 
breathing 0(tt, 0) phonon mode. Cul, Cu2, Cu3 and 
Cu4 label the copper sites, and 01, 02 and 03 differ- 
ent oxygen sites. For the covalent el-ph coupling, the 
hybridization term t l ^ d is modified (a / and f3 = 0), 
whereas for the ionic el-ph coupling, the Cu-site energy 
e l d is changed (a = and (3^0). 

To clarify different effects of lattice displacement, we 
examine the u dependence of 
(i) the electronic energy 
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(ii) the kinetic energy 
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Here, 



means the ground-state expectation value. 



FIG. 2: (a) Change AE (per unit cell) in electronic energy vs 
the lattice displacement u in the hole-doped case; (b) changes 
AE a (per unit cell) in kinetic energy (open circles), O-site 
energy (open up-triangles), and Coulomb energy (open down- 
triangles) vs the lattice displacement it for the covalent el- 
ph coupling in the hole-doped case; (c) same as (b) for the 
ionic el-ph coupling, and open squares represent the Cu-site 
potential energy; (d), (e), and (f) same as (a), (b), and (c), 
respectively, in the electron-doped case. 



In Figs. Ebi and we show the changes of dif- 
ferent energies as a function of the lattice displacement 
u in the hole-doped case. The corresponding results in 
the electron-doped case are shown in Figs. 03, and 
|2f. In both hole- and electron-doped cases, the electronic 
energy is reduced by lattice displacements. As one can 
see, in the hole-doped case, the covalent and ionic el-ph 
couplings have an equal effect in the reduction of elec- 
tronic energy, while in the electron-doped case, the ionic 
el-ph coupling has a stronger effect than the covalent el- 
ph coupling. 

From Figs. [2b and[2t, for the covalent el-ph coupling, 
one can see that the reduction of electronic energy comes 
from the gain in kinetic energy, while for the ionic el-ph 
coupling it originates from the gain in Cu-site potential 
energy as shown in Figs.[2t and[2f. Notice that although 
the noninteracting model (Ud — 0) also decreases its elec- 
tronic energy by lattice displacements through the same 
mechanism, i.e., the gain in kinetic energy or Cu-site po- 
tential energy, electronic correlations play an important 
role in the strongly correlated system. As a matter of 
facts, for a given lattice displacement u the Coulomb in- 
teraction Ud transfers some charges from copper sites to 
oxygen sites, making the kinetic energy Ek in the inter- 
acting system even lower than that in the noninteracting 
system. This demonstrates that electronic correlations 
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FIG. 3: (a) Charge densities at the copper sites Cul (filled up- 
triangles) and Cu3 (filled down-triangles) vs u for the covalent 
el-ph coupling in the hole-doped case; (b) same as (a) but 
at the oxygen sites Ol (open up-triangles) , 02 (open down- 
triangles), and 03 (open circles); (c) and (d) same as (a) and 
(b), respectively, in the electron-doped case. 



FIG. 4: Spin correlations between the nnn Cu sites vs u for 
the covalent el-ph coupling in the hole-doped case; (b) same 
as (a) for the ionic el-ph coupling; (c) and (d) same as (a) 
and (b), respectively, in the electron-doped case. The open 
circles, open up- and down-triangles represent the site pairs 
(i,j)=(Cul, Cu3), (i,j) = (Cul, Cu2), and (i,j) = (Cu3, Cu4) (see 
Fig.0, respectively. 



contribute essentially to the gain in kinetic energy. In 
addition, we observe that for the ionic el-ph coupling 
and in the noninteracting model, the Cu-site potential 
energy decreases more in the hole-doped case than in the 
electron-doped case. On the other hand, Figs.Et and|2t 
show that in the strongly correlated system this behavior 
is reversed, which is consistent with the physical situation 
where the doped electrons mainly go to the copper sites, 
as a result, the charge fluctuations at the copper sites are 
much stronger than those in the hole-doped case. 



have suggested that the quasiparticles of the low-energy 
excitations in an antiferromagnet can be well described 
by the doped hole or electron dressed by local antifer- 
romagnetic spin fluctuations (spin-bag picture). There- 
fore, it is important to understand the effects of lat- 
tice displacements on local spin correlations in order to 
shed light on the single-particle excitation as measured in 
ARPES. Here, we are interested in the spin correlations 
of the next-nearest-neighbor (nnn) Cu-Cu sites, which 
are defined as 



In general, the strong el-ph coupling in both p- and n- 
type cuprates can be attributed to the strong interactions 
of lattice displacements with local metallic charge fluctu- 
ations. McQueeney et al. 4 have studied the doping de- 
pendence of phonon densities of states in ha2- x ST x CUO4 
and found that the development of the 70 meV phonon 
mode near the doping-induced metal-insulator transition 
is due to the low-energy one-dimensional charge dynam- 
ics. Moreover, Ishihara et a/piSiii have studied the one- 
dimensional two-band Peierls-Hubbard model and indi- 
cated that the charge-transfer fluctuations are strongly 
enhanced near the transition from a Mott insulator to an 
ionic insulator (or a metal). As shown in Fig. El due to 
lattice displacements, the charge densities at both cop- 
per and oxygen sites are modified in a similar way for the 
hole- and electron-doped systems. However, in our cases 
the charge-transfer fluctuations have little contribution 
to the covalent el-ph coupling, which is evidenced by a 
rather small change of the O-site energy (see Figs.l^t and 

120- 

Now we turn to discuss the effects of lattice displace- 
ments on the magnetic properties. Previous studies^i^ 



S(M) = (Si-Sj), (7) 
with the spin operator Si = J^ a p AoPaP^W- 



Figure 0] shows the nnn Cu-Cu spin correlations as a 
function of u. One can clearly see that for both hole- and 
electron-doped cases, the el-ph couplings have significant 
influences on the nnn Cu-Cu spin correlations: while the 
spin correlations in the direction parallel to the lattice 
displacement have little change, the spin correlations in 
the direction perpendicular to the lattice displacement 
increase or decrease alternatively depending on both the 
Cu sites and the doping. In particular, there exists a 
crucial difference between the hole- and electron-doped 
cases. In the hole-doped case, the spin correlations in- 
crease for the Cu sites with smaller charge density, and 
decrease for the Cu sites with larger charge density, while 
in the electron-doped case, this behavior is reversed. 
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FIG. 5: Same as Fig. ^]for the normalized spin correlations. 

In order to filter out of the effects of charge fluctua- 
tions, we have divided the spin operator by the probabil- 
ity of single occupancy. The corresponding normalized 
spin correlations are defined as 

S*(U) = f^. (8) 

where P S i =< n% — 2njj • mi > denotes the probability of 
single occupancy at the site i. Eq. ©, thus measures the 
spin correlations normalized to the spin density. Figure|S] 
shows the normalized spin correlations as a function of u. 
From Figs. and[SJ}, it is clear that in the hole-doped 
case, both ionic and covalent el-ph couplings induce a 
strong response of spin degree of freedom to phonons. 
On the other hand, we observe a rather weak effect in 
the electron-doped case, which is shown in Figs. Et and 
Hi. 

The above results indicate that there exist dramatic 
differences of the spin-phonon coupling in p- and n-type 



high-T c cuprates: the coupling of spin with lattice dis- 
placements is much stronger in the p-type materials than 
in the n-type materials. This could give an explanation 
why ARPES 2 detected a strong mass renormalization 
in the p-type cuprates, while it is absent in the n-type 
cuprates. The reason is that within the spin-bag picture, 
the quasiparticle property is mainly determined by local 
antiferromagetic spin fluctuations and a strong change of 
local spin correlations can dramatically affect the motion 
of doped charge carriers in the background of antiferro- 
magnet. We argue that if the similar difference in spin- 
phonon couplings occurs for other phonon modes at small 
and large momenta, we then expect that in the n-type 
superconductors, the phonon contributions to electron- 
electron pairing and electronic transport should be much 
smaller than those in the p-type materials. 

Summarizing, we have presented QMC studies of the 
effects of the half-breathing 0(ir, 0) phonon mode on 
the ground-state properties of the 2D three-band Peierls- 
Hubbard model. We found that in the physically inter- 
esting parameter regime, the el-ph couplings have similar 
effects on the ground-state energies and charge fluctua- 
tions in both hole- and electron-doped systems. This 
suggests that the strong softening of the half-breathing 
O(7r,0) phonon mode in both p- and n-types materials 
can be attributed to the the strong interactions of lattice 
distortions with local metallic charge fluctuations. In ad- 
dition, our results show that the spin-phonon coupling is 
much stronger in the p-type cuprates than in the n-type 
cuprates. This difference in the spin-phonon coupling 
could provide an explanation for the difference of mass 
renormalization observed in ARPES experiments. 
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